Helixing Transition in Filaments under Torque 
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We show that a filament in thermal equilibrium subjected to torque r and a large enough tension 
undergoes a phase transition of infinite order when the torque is twice the temperature T: for very 
large values of the ratio t/T, the polymer's average supercoiling responds linearly to the applied 
torque, but it tends to zero exponentially fast as the ratio reaches the critical value t/T — 1/2. 
When approaching the transition from below, correlation lengths tend exponentially to infinity, 
which suggests an essentially topological nature. At room temperature, this transition happens at a 
torque ~ 2 pNxnm, and it should be detectable in carefully performed single molecule manipulation 
experiments. 

PACS numbers: 82.35.Pq,05.70.Fh, 87.14.E-, 87.14. G- 



The thermodynamics of filaments has been the center 
of extensive study over the last four decades, since the 
pioneering works of Edwards, De Gennes and others [H[2] 
started what is now called soft condensed matter physics. 
Motivated by practical reasons and advances in polymer 
science, by e.g. Kuhn and Flory [3J, research in polymer 
physics initially concentrated on the study of collective 
statistical behavior: solution, gelation, networks, entan- 
gling and unentangling, rheology, etc. [2]- 

Biological and technological interest in highly struc- 
tured polymers, most notably DNA, have motivated 
experimental research into their behavior as single 
molecules under precisely controlled mechanical manip- 
ulation. Direct single molecule manipulation of DNA, 
proteins, and other more or less complex polymers [4] 
has opened the field to the analysis of entirely new phe- 
nomena, such as writhing and load/temperature induced 
denaturation. Meanwhile, advances in nanotechnology 
have shown how single polymers and DNA can be com- 
pose into novel materials [5], whereas novel inorganic 
nano-filaments such as carbon nanotubes were also in- 
troduced [5J. The confluence of theoretical, experimen- 
tal, and technological developments has made it timely to 
analyze the behavior of single (rather than many) poly- 
mers under load; much effort [THTTj has concentrated on 
understanding single molecule experiments for DNA in 
particular [HHUJ- 

Here we show that an infinitely smooth helixing tran- 
sition is expected in a long filament under applied load 
in a thermal bath. The transition is controlled by the 
ratio between applied torque r and temperature T (mea- 
sured as energy): for T <C t, the polymer responds lin- 
early to torque as it assumes an helical structure. For 
T > r the order parameter (helicity) disappears and the 
structural response of the polymer to torque is cancelled 
by thermal fluctuations, as its generalized helical rigidity 
becomes infinite. The transition is of infinite order and 
characterized by an exponential divergence of correlation 
lengths. 

Consider modeling the following experiment: a meso- 



scopic one-dimensional system, which for simplicity we 
take as a polymer with monomers of length a, is held un- 
der tension x an d subjected to a torque r by magnetic or 
optical tweezers [4]. Held straight, under applied torque 
it will respond with an elastic change of the internal tor- 
sional angle. On the other hand, if the torque is large 
enough, one expects a structural change, in which the 
polymer assumes an helical configuration. 

We consider a continuum energy model for this system. 
Suppose the tension is large enough that x - 2 exceeds the 
flexural rigidity |21j : then the dominant contribution to 
the energy will come from the change x^/i in distance h 
between the ends by which the filament is held under ten- 
sion x- If? i n a continuum limit, ip(s) = r(s) exp [i#(s)] 
represents the deviation from the straight filament config- 
uration in the plane perpendicular to it , and s is the in- 
trinsic coordinate (see Fig. 1), then the contribution (per 
unit length) to the energy density from tension, neglect- 
ing subdominant terms, is —xdh/ds — x(r 2 + r 2 9 2 )/2. 
Here dh is the experimentally measured change in dis- 
tance h between beads. The contribution from the torque 
is — rt ■ d, where t is the unit vector tangent to the curve 
and d the Darboux vector. For configurations that wind 
around the experimental axis, it can be approximated as 

-T0. 

One arrives at the following expression for the energy 
per unit length 
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(p 2 + P 2 uj 2 ) - g{p)ru, 



(1) 



where p £ [0, oo), w = 6e (— oo, oo) are functions of the 
parameter s £ [0, L), with L large in the thermodynamic 
limit. The regularization at small p, g(p), is g — 1 when 
p > Po (a cutoff at small p, po is of the order of the lateral 
size of the polymer), and g — otherwise. A less rigorous 
choice does not affect the following conclusions. The cut- 
off g(p) avoids unphysical configurations in which p = 
yet u^O. More importantly, it gives a lower bound to 
the energy II]) . Physically it describes the breakdown of 
the model when the radius of the helix becomes compara- 
ble to the radius of the polymer, and couplings between 
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FIG. 1: The polymer held at its ends by a tension \ along 
the direction h and subject to a torque r assumes a helical 
configuration (top right): the two-dimensional vector tp(s) 
(red arrow) describes the deviation from the straight filament 
configuration in the plane perpendicular to h, while s is the 
intrinsic coordinate. For flexural modes of the helix (bottom 
right, blue arrow), which are suppressed by tension and do 
not couple with the torque, tp describes the deviation from the 
flexural mode. Top left: as t/T — > 1/2+ the average helicity 
(uj)1 (red solid line) tends to zero exponentially fast, while 
both the generalized helical rigidity (k) = l/{k~ 1 } (dashed 
black line) and the correlation length £ (black solid line) ap- 
proach infinity with exponential behavior. Bottom left: in an 
experiment of single molecule manipulation one often mea- 
sures the torque after having fixed the supercoiling a. The 
transition is then revealed by plotting the torque vs. o — r / p 
to suppress the inner torsional variable of rigidity p. 



the helical gain and the inner torsional degree of freedom 
become relevant. 

One can check that, if p is kept fixed, the energy in (fTl) 
is minimized by an helix of gain angle per unit length 



oj = k(p) i 



(2) 



where k(p) is the linear density of helical rigidity of the 
structure per unit length (and therefore k~ l is the helical 
compliance) and normalized to temperature. We have: 
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The energy of the helix is H[p, ui(p)] 

r 2 g( P ) 



V{p), with 



V{ P ) 



2 X p 2 



(4) 



If there is no constraint on the radius p, the low- 
est energy structure corresponds to p = po an d helicity 
ujq = k(po)~ 1 T. More realistic choices for the cutoff g(p) 
accommodate the statics of different polymers but have 
no effect on the transition. 

Let us now investigate the effect of thermal fluctua- 
tions on this structure. The term —tui disregards the 
energy of polymer configurations winding around a dif- 
ferent axis than the experimental one: nonetheless, one 
can always decompose these configurations in a slowly 
varying flexural mode-suppressed by tension-plus a he- 
lical mode described by our hamiltonian. The flexural 
modes (Fig. 1) are suppressed by the tension and do not 
couple to the torque, therefore representing an inconse- 
quential background for our torque-driven transition. 

In the presence of thermal fluctuations, entropy reduc- 
tion is the cost of structure. Helical fluctuations , al- 
though favored by energy, should be entropically disad- 
vantageous compared to planar ones. One can anticipate 
that this competition between entropy and energy might 
drive a transition. 

We can gain some heuristic insight by computing the 
free energy of a helix of fixed angular gain Q. We sum 
over fluctuations in p only while maintaining Q fixed, 
to obtain the partition function of an helix of gain u), 
or Z Q = e-^ J[dp] e -^ti(p 2+p2Q2 ) ds , where [dp] = 
EL=i dp(ai) in the limit of small a. This path integral 
can be computed exactly by turning it into the trace of 
the density operator for the harmonic Schrodinger prob- 
lem in p, which then projects on its lowest eigenvalue 
|o)|/2 in the thermodynamic limit of large L. Thus one 
finds the free energy density f a = L^ 1 hiZ for an helix 
of fixed helicity uj. By subtracting the equipartition free 
energy density (which is Tin (x/T) /2a, obtained by set- 
ting Q = 0) one finds the extra free energy contributed 
by a structure of defined helicity ui, 

Af = -TU + T\Q\/2, (5) 

which corresponds to an energy reduction -tw, but also 
an entropy reduction As = — |cj|/2. Their interplay 
presages a transition: for values of \t\/T < 1/2 helical 
structures are suppressed, as any helicity w/0 gives a 
positive contribution to the free energy. However when 
\t\ > T/2, or equivalents T < T c = 2|r|, the entropy 
cost of helicity can be offset by energy gain and helical 
structures of the same orientation of r lower the free en- 
ergy [22]. As we will see below, this heuristic argument 
correctly predicts the critical temperature T c = 2|r|. 

To obtain the exact nature of the transition we return 
to the full partition function 



(3) 



Zoc f[pdujdp}e- l3 fa' Hds . 



(G) 



Since H in (fTl) is quadratic in uj, the helicity can be inte- 
grated out. This removes the p factor from the functional 
measure in (raj) and provides an equipartition factor, re- 
turning 



Z oc 
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Y" f[dp]e-^o[iP 2 +y(p)}^. (7) 



From (rol) we can define the order parameter for the he- 
licity 



(w) = TL~ l d T In Z= (fe _1 )r, 



(8) 



where (k^ 1 ) is the generalized helical compliance. Equa- 
tion pi) can be rewritten as (k^ 1 ) = 2TL~ 1 d T 2 \nZ and 
then from fff\ and Q one obtains 
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the generalized helical compliance is simply the thermal 
average of the "static" helical compliance given in (pf. 

The helical compliance is non-zero when the radius is 
finite. From p| and Ety, when (w) =^ 0, we can picture 
the system as fluctuating around a helix of gain angle (uj) , 
and of radius given by (p~ 2 ) = x(k~ x ) = (&)/t: consis- 
tent with our heuristic findings, the radius is confined 
when helical thermal fluctuations dominate. 

The partition function in fcfh can be expressed through 
standard path integral manipulation in terms of a trace 



Z \ | - | Tr ! < 
.X 
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of a one dimensional Schrodinger operator 
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(10) 



(11) 



[The factor of proportionality in ( 10 1 depends on bound- 
ary condition and is non-extensive in L. In the absence 
of truncation g the operator in ( |11[ ) is unbounded from 
below, just like the energy in (fTl).] 

In the thermodynamic limit of large L, the trace in 



( 10 1 projects on the lowest bound state of eigenvalue 
eo — if there is one. The pre-factor {T/x) L in (10 1 only 



contributes two equipartition terms corresponding to the 
two oscillating variables u) and p and is independent of 
t. Thus one has (w) = —Td T e ^ 0, giving rise to he- 
lical structures of confined radius, as described above, 
whenever there is a discrete spectrum. 



The operator in (11) corresponds to the well known 



Schrodinger problem of an attractive truncated inverse 
square potential. Its discrete spectrum is known to dis- 
appear when |t|/T < 1/2, which confirms the existence 
of a transition at T c = 2|r|. Indeed since the contribution 



from the continuum spectrum in ( 10 ) is non-extensive in 



L, for T > T c the partition function reduces to equipar- 
tition alone, or Z = (T/xY L ^ a \ which is independent 
of r, and thus (w) = above the transition: the helical 
structure is lost. 

When T — » 2t~ the 1-D Schrodinger problem of a 
truncated inverse square potential on a half line has 
been studied by Gupta and Rajeev [15 through renor- 
malization methods, later extended to general dimension 
by Camblong et al. |16) . Indeed the problem is a long 
standing one of quantum mechanics, and has relevance 
in molecular physics [17J , string theory |18j , and polymer 
physics (in different context than ours) [T5] , 

Defining v 2 = (t/T) 2 — 1/4, the remaining bound 
eigenvalue can be approximated in the limit v — > + |15j 
or t/T -»■ (1/2)" 



as 



e = -I 



f [l + OM], 



(12) 



with I = p\e 7 x/T, and 7 the Euler-Mascheroni con- 



stant. From (12 1 it follows that the helical order param- 
eter 
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(13) 



disappears at transition with infinite smoothness (all its 
derivatives tend to zero), Fig. (1). The transition is 
therefore of infinite order, as one might anticipate given 
its topological nature. Clearly, at transition the general- 
ized rigidity t/(oj) approaches infinity exponentially fast 
and therefore, from (pM), so does (p 2 ). 

Finally, it is interesting to reason in terms of helical 
pitch, 2it/(lj). In particular its ratio with the radius of 

Therefore the pitch of 



the helix, from fl9| is ex (uj 



'1/2 



the helix in units of the radius of the helix approaches 
infinity exponentially fast at transition, i.e. the helix be- 
comes degenerate to a plane infinitely fast. Above the 
transition, as we have seen, only planar fluctuations con- 
tribute to the partition function. 

Our previous heuristic treatment lead us to suspect 
that fluctuating helical configurations provide a domi- 
nant contribution to the partition function below the 
transition. In contrast, only planar oscillations con- 



tribute above the transition. Eq. ( 13 ) shows that these 



helical structures disappear very smoothly. 



From ( 10 ) one has for the free energy density / = 



Te Q + a~^Tln(x/T). Neglecting the last term, which 
corresponds to equipartition, one finds from s = —dxf 



As 



1 



(uj) [1 + 2z/ 2 + 0(u 3 )} 



(14) 



which shows that close to transition, the principal part 
of the entropy reduction is proportional to — (w)/2. This 
is the same result we obtained from ^ for the struc- 
tural cost in entropy for a helix of uniform gain, sug- 
gesting that at transition the order parameter tends to 
uniformity as it disappears: its fluctuations tend to zero, 



while their correlation length approaches infinity. The 
first statement is obtained immediately by differentiat- 
ing the expression for (ui) in ( f 3 ) with respect to r. We 



demonstrate the second one below. 

The infinite order character of the transition as well 
as its topological nature should also be evidenced in an 
exponentially smooth disappearance of the correlation 
lengths. These can be computed as follows. We intro- 
duce a varying external field t(s) = r + r](s) with r 
uniform, as before. Then we compute the functional 
Z[ri(s)] by summing over all the trajectories in oj(s), 
p(s) as in J6J. Clearly the helical correlation function 
G(s 1 ,s 2 ) = (w(si)w(s 2 )) ~ ( w ) 2 is given by [H] 



G(si,s 2 )=T 2 



8 2 In Z [rj\ 
Srj(si)5r](s2) 
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(15) 



The new partition function Z[rj\ for rj(s) ^ is still 



given by a trace as in (10 1, but with the replacement 

LH -> LH + 4> fa r)(s)ds, with cf>(j>) = gT/(T X p 2 )- Us- 
ing imaginary time s— dependent perturbation formal- 
ism, one obtains through standard calculations 



rp2 roo 

G(si )S2 ) = — / dfc|<c |^|fc>] 



2 e -Ae fci0 |si-s 2 | 



(16) 



where \k) is the unbound eigenvector of ( 11 1 of eigenvalue 
Tk 2 /(2 X ) and Ae kfi = Tk 2 /(2 X ) - e 0) where e is the 
remaining bound eigenvalue. 



From (16), by bringing exp(eo|As|) out of the inte- 



gral and noticing that exp [— T\As\k 2 (2x) 1 ] > 

<5(/c)(27rx/3/|As|) 1 / 2 , for large |As|, one finally has 

O eo|si-s 2 | 



G(s 1 ,s 2 ) oc 



\/ s l - s 2 ' 



(17) 



showing that the correlation length £ = — 1/eo ap- 
proaches infinity exponentially fast at transition: 



£~ exp 
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V(t/T)2-1/ 4/ 



(18) 



Some similarities are evident in this behavior with the 
Kosterlitz-Thouless (KT) transition, but this should not 
be over interpreted, as we are dealing with an externally 
driven situation. Unlike the KT case, here we do have 
an order parameter, (w), which disappears at transition. 
Also, the external field, r, breaks the symmetry of the 
problem. In fact, if we want to establish a formal sim- 
ilarity between the two transitions, we must note that 
in the formula for the correlation length r plays the role 
of T in the KT transition. Similarly, while our heuristic 
argument in ([5]) is reminiscent of the analogous one for 
KT, energy and entropy are exchanged: while in KT the 
formation of an unbound topological structure, the vor- 
tex, is favored by entropy at high T and disfavored by 



energy, the formation of a helical structure is favored by 
energy at high r and disfavored by entropy. 

In experiments of single molecule manipulation, one 
fixes the supercoiling and measures the resulting torque. 
The transition proposed here can thus be revealed ex- 
perimentally by reducing the supercoiling to zero and 
observing that the torque converges very slowly to the 
value T/2, which at room temperature corresponds to 
^2 pNxnm. In practice there will also be an inner tor- 
sional degree of freedom £ of stiffness /i. Therefore the 
total supercoiling will be given by a = (u>) + r/ fx: exper- 
imentally the transition will be manifested by plotting 
t/T vs. a — r/fi, as shown in Fig. 1. 

In conclusion, inspired by a single molecule problem, 
we have introduced a simple model of a free hamiltonian 
whose helicity is coupled to a uniform external field, and 
shown that it leads to an infinite order topological tran- 
sition. This transition is driven by the competition of 
entropy and energy, both of which are lowered by helical 
fluctuations. 

Our formalism is more general and might be ap- 
plied, for instance, to kinetic transitions in problems of 
Langevin dynamics in an external field w. Indeed a 
path integral formalism applies to that case, with la- 



grangian (x — w) 2 /AD — divw/2 



One can verify 



that the choice of a vortex w — rez A x/|x| 2 + x<j){\x\), 
where r is now the vorticity, can lead to an expression 
similar to the one in M when x is represented in po- 
lar coordinates, and for appropriate choice of the field 
4>. That would lead to a kinetic transition in t/D, due 
to a competition between vorticity and diffusion. Also, 
% can be rewritten in terms of a boson condensate as 
(3H = x\iJA 2 / 2 + iT 9 U>*i> ~ V^) / 2 \^P\ 2 - We will dis- 
cuss these different contents elsewhere. 
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